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Abstract. This article reviews the use of differential forms and Lie derivatives to find sym- 
metries of differential equations, as originally presented in Harrison and Estabrook 1971 PP. 
, An outline of the method is given, followed by examples and references to recent papers 

^ ' using the method. 

o ; 

_ Key words: symmetries; differential equations; differential forms 

00 ■ 2000 Mathematics Subject Classification: 22E70; 34A26; 35A25; 35A30; 53Z05 

I ! 

^ '. 1 Introduction 

^ : 

' In 1969-70, Frank Estabrook and the present author found a method for finding symmetries 

, of differential equations using differential forms and Cartan's formulation of differential equa- 

l/^ I tions jj. (This will be called Paper I.) It was not something we were searching for; rather, we 

■ were simply trying to understand how the symmetries of Maxwell's equations could be found 
. from the differential form version of those equations. Once we realized that the key to sym- 

metrics was the use of the Lie derivative, it became clear how to apply this to all differential 

■ equations. An outline of the method, with examples, will be given here. A few computer prog- 
. rams that use it will be mentioned, along with a number of published papers on symmetries 

^ I which have used it. 

^ ' The reader may wonder at the order of the names of the authors on our original paper. Since 

we had done roughly equal amounts of work on the research, the order was determined by the 
flip of a coin. 

\ The method proceeds as follows. We consider a set of partial differential equations, defined 

on a differentiable manifold M of n independent variables and m dependent variables. (Ordinary 
differential equations constitute a special case; we will mention those later.) We define the partial 
derivatives of the dependent variables as new variables (prolongation) in sufficient number to 
write the equations as a set of first order equations, thus extending the manifold to a manifold M' . 
Then we can formulate those as differential forms. We speak of the set of forms, representing 
the equations, as an ideal /. It is to be closed. 

We can recapture the original set of equations by two operations: we specialize to a submani- 
fold by letting the dependent variables be functions of the independent variables ("sectioning") 
and then we set the pullback of the forms to zero ("annulling"). The resultant equations are 
the original set of first order partial differential equations, and the submanifold is a solution 
submanifold. (All calculations are local and we do not use any manifold structure except diffe- 
rentiability.) 

Lie derivatives of geometrical objects, like tensors, are associated with symmetries of those 
objects. The Lie derivative of a geometrical object carries it along a path, determined by 
a vector v, in its manifold. If the Lie derivative vanishes, then the vector v represents the 
direction of an infinitesimal symmetry transformation in the manifold. A differential form is 
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a type of tensor (totally antisymmetric on the indices of the components), so it has a Lie 
derivative. We may construct the Lie derivative (symbolized by £) of the forms in the ideal I. 
Setting the Lie derivative of these forms equal to zero should therefore represent symmetries — 
except for one thing. When we make an infinitesimal transformation away from the original 
variables, we require that the new form of the differential equations should vanish — but the 
old form must also vanish. Thus we want the original forms in / to vanish, but also their Lie 
derivative must vanish when that happens. In other words, we require that the Lie derivatives 
of the forms in / to be linear combinations of those forms themselves — and when they vanish, 
then the Lie derivatives also vanish. We can express this by writing = (mod I), or 

C /. (1) 

This is satisfied by letting the Lie derivative of each differential form in J be a linear combination 
of the forms in /. 

Equation (pQ) will contain a number of Lagrange multipliers, the coefficients of the forms in 
their linear combinations. Those are to be eliminated. Once they are eliminated, there remains 
a set of linear homogeneous first order equations for the components of v in M', which are the 
symmetry generators. The equations are simply the determining equations for the symmetries 
of the original set of differential equations, considered as point transformations in M' . 

We note, in this set, that the derivatives of the components of v (the generators) will be taken 
with respect to both dependent variables (including prolonged ones) and independent variables. 
One often assumes that the generators for the independent variables (often denoted by ^ and rf) 
are functions only of the independent variables (there are exceptions.) The determining equa- 
tions will usually show that feature promptly. 

Some examples of familiar equations will be presented to show how the method works. Some 
of this material was presented by the author at the second Kiev symmetry conference in 1997 
and can be found in its Proceedings [2j. This will be denoted as paper II. 

In paper I, we adopted the term "isovector" for the vector v, describing a symmetry trans- 
formation in M', even though the term had been used elsewhere in the physics literature. We 
did not think there would be any confusion. A number of authors thus refer to this method as 
the "isovector" method. 

2 Lie derivatives of differential forms 

First we note some simple features of Lie derivatives of differential forms. (See paper I.) 

(1) Lie differentiation preserves the rank of a form. 

(2) The Lie derivative of a coordinate is simply the component of v in that direction: 

(3) The Lie derivative of a function on M' (0-form) is simply its directional derivative: 

= v(/) = ^;7,^• 

(Commas represent partial derivatives. Sometimes they will omitted when the context is 
clear.) 

(4) The Lie derivative of a wedge product obeys the Leibniz rule (the subscript v may be 
suppressed where it is not necessary): 



£{a Af3) = {£a) A /3 + a A {£I3). 
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(5) The exterior derivative d and the Lie derivative £ commute. In particular, 

£vdx^ = d{£vX^) = dv^; 

the Lie derivative of the differential of a variable equals the differential of the corresponding 
component of v. 

3 The one-dimensional heat equation 

We write the one dimensional heat equation 

Uxx = ut (2) 
as a first order set of equations by defining a new variable w: 

Ux =w, Wx = Ut- (3) 
The variables are x, t, u, and w. We construct two 2-forms by inspection: 

a = du Adt — wdx A dt, 

(3 = dw A dt + du A dx (4) 

(a is a contact form.) If we "section" these forms — specialize to a submanifold u = u{x, t) and 
w = w^x, t) — we get 

a = {uxdx + utdt) Adt — wdx Adt = {ux — w)dx A dt 

and 

P = {wxdx + Wtdt) Adt + {uxdx + Utdt) A dx = {wx — Utjdx A dt 

where we have used the antisymmetry of 1-forms. We now "annul" these forms — set them 
equal to zero — obtaining Eqs. the original first-order set of equations. The forms a and /3 
in Eqs. ^ now constitute the ideal / of forms representing the heat equation @. 

We note that / is not unique; we may as well represent the heat equation by defining z = ut 
and constructing an ideal /' with a 1-form 

7 = —du + wdx + zdt, 
its exterior derivative 

d'y = dw A dx + dz A dt, 
and the 2-form 

5 = dw Adt — zdx A dt. 

We note that a = —7 A dt and /? = 5 — 7 A dx. 

We work first in the ideal /. Write the Lie derivatives of a and /3 as linear combinations of 
themselves. Expand the Lie derivatives by the rules above. We also drop the wedge product A 
and the subscript v on to save writing. 

£a = £{dudt — wdxdt) 

= {£du)dt + du{£dt) — {£w)dx dt — w{£dx)dt — wdx{£dt) 
= dv^dt + du dv^ - v^dx dt - wdv^dt - wdx dv^ 
= \i{dudt — wdx dt) + A2 {dw dt + dudx). 
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The Aj are 0- forms (functions). Expand the dv^ by the usual chain rule, since the are functions 
in M' , using all four variables. Since dtdt = 0, etc., by the antisymmetry of 1-forms, some terms 
drop out. We have 

(f ^dti + v^^dx + v'^^dw) dt + du {v^^^dt + v^^^dx + v^^^dw) 

— v^dx dt — w{v^^dx + v^udu + v^^jdw) dt 

— wdx {v\dt + v^^du + v^^dw) 

= \i{dudt — wdx dt) + A2 {dw dt + dudx). 

There will be 4!/2!2! = 6 basis 2-forms {dxdt, dxdu, dxdw, dtdu, dtdw, and dudw.) We equate 
the coefficients of these 2-forms to get 

v^-v"" -w{v% + v\t) = -wXi, 
-v% - wv\^ = -A2, 
-wv^w = 0> 

-v'^u-v\t + WV% = 

v% = 0. 

Eliminating the Lagrange multipliers Aj gives us one half of the determining equations: 

Expansion ol £(3 gives us the other half. One quickly sees from them that is a function of t 
only and that is a function only of x and t. Further calculation gives the usual six generators 
plus addition of an arbitrary solution. Exponentiation of the transformation proceeds by setting 
V • 7 = (contraction of v and 7, symbolized by a dot) and solving, where w and z are replaced 
by their values as derivatives of u (the usual method) . 

Another way to proceed, which removes the need for the multipliers, is to use a = to 
replace dudt, anywhere that combination occurs in the expansion of the Lie derivatives of a 
and /3, by wdxdt, and to use /? = to replace dw dt by —dudx. This may save considerable 
work in complicated cases, especially in cases where not all forms in the ideal are of the same 
rank. In those cases, some of the Lagrange multipliers may need to be forms (of rank greater 
than zero) themselves in order for the right hand sides to be of the same rank as the left hand 
sides, and that means that there may be very many coefficients to be eliminated. If one can 
avoid that, labor may be saved. 

One can also use the ideal /' for the heat equation. (This is the technique used in paper I.) 
There are now five variables in M': x, t, u, w, and z. In I', there is only one 1-form 7, and so 
its Lie derivative equation is simple: 

£1 = A7, (5) 

where A is a multiplier. One can expand the Lie derivative by an identity for any form uj, using 
the contraction operator: 



£v(jJ = d(y • (j) + V • duj. 



(6) 
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Write F = V • 7, which is a function, and expand Eq. © using Eq. @ (with to = j ) and 
identities for contraction (e.g., v • {dxdy) = v^dy — v^dx) (see paper I). We get 

F = + wv"" + zv^ (7) 

and 

= dF + v'"dx - v^dw + v'^dt - v^dz 
= A7 = X{—du + wdx + zdt). 

Expand dF with the chain rule, equate coefficients, eliminate A and use Eq. Q, and we get all 
generators f* in terms of F and its derivatives (subscripts on F are derivatives): 

v"" = F^, v' = F,, v'^ = -F + wF^ + zF,, 
v"" = -F^ - wFu, v' = -Ft - zFu. 

The exterior derivative of Eq. © is 

d{£-i) = £{d-f) =dXA-f + Xd-f 
so that the Lie derivative of dj is also in the ideal. There is now only one equation left: 

£S = Xi6 + A2d7 + T A 7, 

where the Aj are 0- forms and r is an arbitrary 1-form with five terms. The term in du in r 
can be eliminated by substituting from 7, and that drops out. This procedure also gives the 
standard determining equations. 

4 Computer programs 

There are a few computer programs which use this technique. Some of these were written, 
in REDUCE, by D.G.B. Edelen The programs are probably still available from Lehigh 
University. Other programs were written by Gragert, Kersten, and Martini, also in REDUCE. 
They published several works which developed and used this software, including a program for 
symbolic integration of overdetermined systems [IJ El El El E] • Problems treated in references 
[TUl im [T^ are studies of a nonlinear diffusion equation 

A(nP+^) + ku" = ut, 

(where A is the Laplacian), the massive Thirring model, and the Federbush model. The present 
author used one of the programs with E.D. Fackerell to explore a relativity problem, in unpub- 
lished work, and one of Fackerell's students, Ben Langton, used it for his Ph.D. dissertation on 
certain solutions of the Einstein equations ^3] . 

Another useful computer program is liesymm, a program found in MAPLE, based on a paper 
by Carminati et al. It works quite well, and there is an additional program called autosimp 
in MAPLE which does some integration of the determining equations, although the integration 
may not be complete. These programs are discussed briefiy in Refs. jl51 116j . A student of the 
author's, David Neilsen, did a master's thesis with liesymm on Einstein's equations |17j . 

An extensive review of symbolic software was done by Hereman |18j in 1997. A nice table 
of programs is provided. No specific distinction is made in that paper between the traditional 
method and the differential form method. 
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5 Nonlinear Boltzmann equation 

This example is actually listed as an example in liesymm in Maple 9.5, although it is not 
worked out. The equation is: 

Uxt + + = 0. 

Possible ideals are /' with five variables, x, t, u, p = Ux, q = ut. 

a = du — pdx — qdt, 

da = —dpdx — dqdt, 

(3 = —dt dq+ {p + u^)dx dt 

or / with four variables, x, t, u, p: 

J = dudt — pdx dt, 

6 = dpdx + {p + v?')dt dx. 

The calculation is quite similar to that for the heat equation. There are four generators. 

6 Vacuum Maxwell equations 

From paper I we write the usual 3-forms that represent the vacuum Maxwell equations in 
rectangular coordinates. Subscripts represent components. 

a = dEx dx dt + dEy dy dt + dE^ dz dt + dB^ dy dz + dBy dz dx + dBz dx dy, 
(3 = dBx dx dt + dBy dy dt + dB^ dz dt — dE^ dy dz — dEy dz dt — dE^ dx dy. 

We simplify these forms by defining ^ = a + if3 and {A, B, C) = (cyclic E^^ + iBjS) = h. Then 
(paper I) we can write 

7 = dh • {dvdt - {l/2)idY x dr) 

or 

7 = dA [dx dt — idy dz) + dB [dy dt — idz dx) + dC [dz dt — idx dy). (8) 

The forms in the ideal will then be 7 and 7*, where the star represents complex conjugate. The 
variables will be t, x, y, z. A, B, C, A*, B*, C*. The generators for the coordinates t, x, y, z 
will be real. The equations for the Lie derivatives are then 

£7 = A7 + ^7* 

and its complex conjugate. 

In paper I the determining equations were worked out by using a vector-dyadic formalism. 
Here we use Eq. (jH| for 7, which is a little clearer. We work with 3-forms in ten variables, 
so that there are 10!/7!3! = 120 different basis 3-forms. There are two equations, for £7 and 
£7*, so that we apparently have 240 equations. However, 120 of them are simply the complex 
conjugates of the others. So we just look at the £7 equation. We see immediately by inspection 
that 3-forms with all terms being (i(field variable) do not appear. Terms of the form (i(field) A 
(i(field*) A d(coordinate) yield only equations for the derivatives of the coordinate generators 
with respect to the complex conjugate fields (which are zero). Thus the coordinate generators 
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do not depend on the complex conjugate fields, nor (by complex conjugation) on the fields. Thus 
they depend only on the coordinates themselves. It is also easy to show that the field generators 
depend only on the fields and not on their complex conjugates. 

This reduces the number of equations to 22-18 with one field 1-form and two coordinate 
forms and four with three coordinate forms. From the first set we get the conformal Killing 
equations plus some expressions for the derivatives of the field generators. From the second we 
get equations for the field generators that are the Maxwell equations themselves. Solution of the 
determining equations gives the 17-generator set given in paper I (15 conformal Killing vectors, 
a scale change on the fields, and a duality change on the fields, plus the addition of an arbitrary 
solution). Steeb presents other symmetries besides these, which depend on the derivatives 
of the fields. (Steeb and collaborators also treat various versions of the Dirac equation |2()j.) 

7 Nonlinear Poisson equation 

We consider the equation: 

Uxx + Uyy + U^z = f{u), 

where f{u) is an undetermined function. Subscripts represent derivatives. We define r = Ux, 
s = Uy, t = Uz- Then 

rx + Sy + tz = f{u). 

The ideal / consists of these forms: 

a = —du + rdx + sdy + tdz, 
da = dr dx + ds dy + dt dz, 

(3 = dr dy dz + ds dz dx + dt dx dy — f{u) dx dy dz. 

There are seven variables. 

We may approach this problem by defining a function H = v ■ a, as we did with the heat 
equation. The Lie derivative of a gives all the generators in terms of H and its derivatives, as 
before. Then the only equation we have left is that for £j3, a 3-form. But equating it to a linear 
combination of q, da, and /? is messy. The multiplier of da, a 2-form, must itself be a 1-form, 
which will have six coefficients (we do not include a term in du, because that can be replaced 
by a, and da A a can be included in the a term.) The multiplier of a must be a 2-form — 
and again we can eliminate du terms because they can be replaced by a, and a A a = 0. But 
that still leaves 15 coefficients. The multiplier of P will be a single coefficient. That totals 22 
coefficients that must be eliminated. 

So we consider an easier way. We define a new ideal /', made up of four 3-forms: /?, adydz = 
{—du + rdx) dy dz, adzdx, and adxdy. The latter three forms are equivalent to a alone. 
The Lie derivative of each 3-form must be a linear combination of all four, thus giving four 
multipliers in each equation to be eliminated. The equations are much simpler; it is easy to 
eliminate four multiplier coefficients in each equation than 22, even though there are now still 
4 X 4 = 16 multipliers. The easiest procedure is to write out the equation for the Lie derivative 
of adydz, eliminate the multipliers to get a set of determining equations and then to permute 
x, y, z (and s, t) cyclically. One quickly gets the result that the generators for x, y, z, and u are 
functions only of x, y, z, u, and the generators for r, s, and t are given in terms of a function 
which is precisely the H defined above, H = v ■ a = —v^ + rv^ + sv^ + tv^. 

The equation for £p now has four multipliers, which are easily eliminated. We find quickly 
that v^, yy , depend only on x, y and z and that they obey the Killing equations. The 
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generators for r, s, and t are written out easily, and one ends up with a single equation involving 
f{u) and f'{u). Solution of that equation for the given / then leads to the final result. 

There is a small technical point. We did not include da in even though we did represent a 
as three 3-forms. Should we have done so? The answer is no. If we take the exterior derivative 
of the 3-forms adydz, etc. we get terms like drdxdydz — in other words, just dadydz. We 
get three of those equations, which are equivalent to the equation for £da. The determining 
equations for / and /' give the same result. 

A similar treatment is used by Satir who writes a set of two-dimensional bosonic mem- 
brane equations as eight 3-forms. He then uses a REDUCE program together with the EXCALC 
differential geometry package to find a 12 parameter group. He remarks that the use of diffe- 
rential forms enabled the calculation to go more quickly that the conventional method. 

8 Nonlinear diffusion equation 

This an equation treated in paper H, originally due to Fushchych — a nonlinear diffusion equa- 
tion with an additional condition. We can write a 1-form as was done above, its exterior 
derivative, and a 4-form for the main field equation. In paper II, it was assumed a priori that 
the generators for the coordinates depend only on the coordinates and that those for the deriva- 
tives of the field were linear in those derivatives. It then turns out that much of the analysis of 
the Lie derivative of the 4-form can be done by inspection. 

9 One-dimensional compressible fluid dynamics 

The equations considered here are (see paper I): 

Pt + ipu)x = 0, 

put + puux + C^Px = 0, 

where isentropic flow is considered so that the pressure is only a function of the density, P = 
P{p), and = dP/dp. The generators for x and t include a p- and u-dependent case, which 
turns out to give the hodograph transformation. 

10 Nonclassical symmetries 

One can generalize the ideal / by including contractions of v with some of the differential forms. 
An example of this was provided in paper I for the heat equation, in which it was shown that 
one gets the equations for "nonclassical symmetries" of that equation, the same equations found 
by Bluman and Cole in 1969 |22]. While this technique has not been explored in detail by this 
author, Webb has studied this set of equations — referred to as a coupled nonlinear Burgers-heat 
equations system — with differential forms and has searched for Backlund transformations for 
the set j2n]- 

11 Ordinary differential equations 

We consider an example: 



y" = fix,y,y'), 
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where the prime indicates differentiation with respect to x. We put z = y' and write two 1-forms, 
a = dy — zdx and (3 = dz — f{x, y, z)dx. There are three variables. The Lie derivative equation 
for a is 

£a = dv^ — v^dx — zdv^ = Xi{dy — zdx) + \2{dz — fdx). 

There are three equations, for the coefficients of dx, dy, and dz; elimination of the multipliers 
gives a single equation, which is an expression for . The Lie derivative equation for /? gives 
another single equation. We assume that the generators and (usually written as and r], 
respectively) are functions only of x and y. In that case, becomes the usual extended generator 
for z = y' and the remaining equation is the usual determining equation for ^ and rj, as given, 
e.g., in Stephani [2l] . 

12 Advantages of using differential forms 

These have been treated in paper II, but are reviewed here. The method is easy to apply. 
One simply writes all equations as first order equations; the differential forms can be written 
by inspection. Calculations may be long because of the necessity of introducing the Lagrange 
multipliers; however, one can choose the ideal to minimize this. One can make use of symmetries 
of form (e.g., cyclic symmetry of coordinates), or one can use the forms to substitute for certain 
terms in the Lie derivative expansion, thus removing the need for multipliers. Independent 
variable generators may easily be considered as functions of the independent variables only (just 
assume that and that simplifies the expansion of the differentials of those generators). 

13 Other examples 

We mention here some research papers in which the differential form method is used. Papachris- 
tou generalized the method to vector-valued or Lie algebra-valued differential forms and treated 
the two-dimensional Dirac equation and the Yang-Mills free-field equations in Minkowski space- 
time |25j (as part of a Ph.D. dissertation with the author.) This was later used to investigate 
self-dual Yang-Mills equations, which work showed connections between symmetry and integra- 
bility (in the form of Backlund transformations) of those equations |261l27l05] . Waller, in three 
similar papers, treats nonlinear diffusion equations (or reaction-diffusion equations) arising in 
plasma physics [25113011^. He uses the technique of writing a 1-form and contracting it with v, 
as done in the second treatment of the heat equation above. 

Edelen has developed the theory of the differential form method extensively. His computer 
programs have already been mentioned. At least two books [HSl |2H] and several papers [HH ESI 
1361 137j explore the use in differential forms in physics, including the method discussed here. 
In papers |35[ I36j he considers a method of characteristics in any number of dimensions, using 
isovector treatments. With this he can write parametric solutions of differential equations. One 
equation he considers is [SHI 

ut Ux = 4ti. 

He gives a solution for the equation as an initial value problem: if u{x, 0) = a{x), with a'{x) ^ 0, 
then 

X = z + a(z)(exp(4r) — l)/a'{z), 
i = (l/4)a'(z)(exp(4r)-l), 
u = a{z) exp(8r). 
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where z is a parameter and r is an arbitrary function. Instead, if one defines r = Ux, one can 
write a simple 1-form for the equation: 

a = —du + rdx + {Au/r)dt. 

Then £a = \a gives equations which yield most of the generators in terms of a function F = v-a, 
which satisfies the linear first order equation (subscripts are derivatives) 

Ft + AFr + (4u/r2)F^ + {8u/r)Fu = AF/r. 

The special solution F = u/r seems to give Edelen's solution, although not all details are worked 
out yet and it is a little uncertain. In Ref. |SZ] he considers "inverse" isovector methods. 

Webb et al consider nonlinear Schrodinger equations for a type of MHD waves, using the 
differential form method |3H]- He also analyzes a nonlinear magnetic potential equation, with 
conservation laws, with the Liouville equation as a special case Pakdemirli and others treat 
boundary layer equations for non-Newtonian fluids, including arbitrary shear stress, power law 
fluids, and other models |40l I41j . §uhubi and others, in a number of papers, consider general 
approaches to equations of balance and other equations EHl El UHl HH HZl HE]- A number 
of these discuss equivalence groups, as a generalization of symmetry groups. One paper with 
Ozer treats nonvacuum Maxwell equations with nonlinear constitutive relations. Another 
discusses steady boundary layer flow past a semi-infinite flat plate j44j . 

Bhutani and Bhattacharya study re-dimensional Klein-Gordon and Liouville equations with 
an interesting approach [JHl- Various types of diffusion equations are treated in Refs. [SOI EH 
1^ . Viscoelastic-viscoplastic rods are studied in Ref. [HSl and power law creep in Ref. 
Equations of meteorology, here meaning steady two-dimensional incompressible inviscid flow 
with a Coriolis term, are studied in Ref. [53 • Hu considers the principal chiral model [SHI ■, using 
differential forms and ideas from Ref. |22j. An interesting paper is that by Barco, who shows for 
a second-order hyperbolic or parabolic differential equation, with one dependent variable and 
two independent variables, that an isovector can be used to generate a similarity solution by 
using a particular Cauchy characteristic vector field j57j . 

Nonlinear thermoelasticity was treated by Kalpakides [SHI- His work is related to that of 
§uhubi [El UHl ISl HH HEl HZl HE] • Hamad and Winternitz considered a generalized nonlinear 
Schrodinger equation, 

izt + Zxx = f{z,z*), 
with attention to both symmetries and Backlund transformations j59j . 
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